Relativity and the lead-acid battery 
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The energies of the solid reactants in the lead-acid battery are calculated ab initio using two different basis 
sets at non-relativistic, scalar relativistic, and fully relativistic levels, and using several exchange-correlation 
potentials. The average calculated standard voltage is 2.13 V, compared with the experimental value of 2.11 V. 
All calculations agree in that 1.7-1.8 V of this standard voltage arise from relativistic effects, mainly from PbC>2 
but also from PbSCU. 
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The lead battery is an essential part of cars, and has numer- 
ous other applications. This well-known invention is now 150 
years old [1 , 2j. About 75% of the World lead production and 
a turnover of about 30 billion USD are due to these batter- 
ies. Although there are electrochemical simulations starting 
from the given thermodynamical data(3j|4]|, we are not aware 
of any ab initio ones for the lead battery. This is in stark 
contrast to other rechargeable batteries, such as the modern 
lithium-ion based systems, where they abound. The problem 
is difficult enough to be a theoretical challenge, and there is 
the additional fascination that, Pb being a heavy-element, rela- 
tivistic effects on its compounds could play an important role, 
as qualitatively found a long time ago|5 - 8 ]. For metallic lead, 
seeEHffl. 

The discharge reaction of the lead-acid cell is 

Pb(s) + Pb0 2 (s) + 2H 2 S0 4 (ag) (1) 
2PbS0 4 (s) + 2H 2 0(/), A£(l). 

The electronic structures of both PbO lfT2ifT4l and f3- 
Pb09 fBlfTTl have been theoretically studied earlier. Es- 
pecially, the metallic conductivity of the /3-Pb0 2 , making 
the large currents possible, was shown to be an impurity, 
conduction-band effect, attributed to donor impurities at oxy- 
gen sites lfT6HT9l . The alloying of the Pb electrode is also 
important in practice, but is not discussed here, because the 
minute amounts of other elements do not affect the EMF of 
the cell. 

The construction of the lead- acid battery! 20] has a posi- 
tive lead dioxide electrode, a negative electrode of metallic 
lead, and a sulfuric acid electrolyte. The discharge reaction 
between a Pb(IV) and a Pb(0) produces 2 Pb(II), in form of 
solid PbSC>4. The experimental thermodynamics of the reac- 
tion are well-known I2T1 . 

The three solids can be treated with existing solid-state the- 
ories, such as density functional theory (DFT), because the 
bonding mechanism in the investigated species is dominated 
by covalent interactions where DFT is expected to provide re- 
liable results. Adequately simulating the liquid phase in multi- 
ple relativistic regimes is beyond current state of the art, how- 



ever. We avoid this problem by introducing the known energy 
AEQ for the experimental reaction 

H 2 0(Z) + S0 3 (#) H 2 SO 4 (0 + A£(§. (2) 

We can use this empirical relationship because only light ele- 
ments and only S(VI) occur in eq. ([2]), whose contribution to 
relativistic effects are small. Combining the equations (1) and 
^ gives 

Pb(s)+Pb0 2 (s)+2S0 3 (#) 2Pb 2 S0 4 (5)+A^(§. (3) 

The voltages for the lead-acid battery reaction may then be 
calculated from the reaction energies 

AE(1) = AE^ - 2AE^ (4) 

where we use calculated AEQ values and experimental 
AE([2]) values. Concentrated sulfuric acid is used in re- 
actions (1) and ([2]). The cell voltage at typical 5.5M (in 
H 2 SO4-10H 2 O) is calculated using the values tabulated by 
Duisman[21 ]. 

Prediction of formation energies in a quantitative manner 
from ab initio calculation requires, in additional to having an 
accurate underlying theory, also absolute convergence of all 
technical parameters and a sufficiently general basis set. It 
is therefore fruitful to approach the problem with several in- 
dependent methods, and see if they converge on the result. 
In our case, we used a linear combination of local orbitals 
(LCAO), with and without a frozen core approximation, us- 
ing the BAND program! 22]; and a full-potential local-orbital 
minimum-basis approach |23l with the FPLO program (ver- 
sion 7.00-27). 

The calculations are performed with crystal structures from 
experimental room-temperature measurements (for structures 
see the Supplementary Material), allowing no ionic relax- 
ations, thus capturing the dynamic electronic effects of rel- 
ativity (meaning Dirac vs. Schrodinger). The alternative, 
more laborious choice would have been to also consider the 
relativistic structural and vibrational changes, including zero- 
point energies (ZPE). The thermal effects on the reaction en- 
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ergies are fairly small and are simply neglected. If a battery 
freezes at low temperatures, it is due to the kinetics, not due to 
the AG. As stated, we assume for simplicity pure, unalloyed 
Pb and a pure /3-Pb0 2 phase in the electrode materials. The 
measured cell voltages for a and /3-Pb0 2 differ by ca. 0.01 
Vl24l. 

In the LCAO approach, the one-electron basis sets repre- 
senting the optimized electron density consist of linear com- 
bination of Herman- Skillman numerical atomic orbitals (NO) 
and Slater- type orbitals (STO). We apply a quadruple-zeta ba- 
sis sets augmented with four polarization functions (QZ4P), 
available in the BAND basis- set repository. Frozen-core ap- 
proximation is consistently applied to reduce the size of the 
variational basis set. The use of frozen core, as implemented 
in BAND[22], is preferable over pseudopotential techniques 
because it essentially allows for all-electron calculations. The 
frozen core orbitals are taken from high-accuracy calculations 
with extensive STO basis sets. For oxygen and sulphur, we 
use all-electron basis sets. For lead we include up to 4f or- 
bitals in the core. The slight change in the deep-core orbitals 
due to formation of chemical bonds was confirmed to be neg- 
ligible. 

The band nature of the solid-state species was inferred by 
studying the density of states (DOS). The orbital character of 
the total DOS was determined with respect to contributions 
from individual atoms. The band structure along a series of 
lines of high symmetry was also analyzed. The overall band 
structures change only slightly with the method: the order- 
ing, the orbital character and the dispersion remain compara- 
ble within LDA and GGA approximations. 

In the case of BAND calculations, we used experimentally 
determined structures both for solid- and gas-phase species. 
The lattice constants and atomic positions were kept fixed in 
all calculations, motived by our interest in purely the elec- 
tronic rearrangement effects due to relativity. The relativis- 
tic effects are investigated by means of the zeroth-order reg- 
ular approximation (ZORA). We consider three cases: non- 
relativistic (NR), scalar relativistic (SR) and the fully relativis- 
tic (FR) with first-order spin-orbit effects taken into account. 
To ensure high accuracy results, the convergence of the cal- 
culations was checked with respect to all crucial numerical 
parameters including the number of k-points, basis set quality 
and size of the frozen core. The formation energies were cal- 
culated with respect to spherically symmetric spin-restricted 
atoms and converged within 1 kJ/mol or less. To sample the 
first Brillouin zone, and evaluate the k-space integrals, we 
used a quadratic numerical integration scheme with 9x9x9 
(Pb and Sn), 5x5x5 (PbO, Pb0 2 , SnO, Sn0 2 ) and 3x3x3 
(PbS04, SnS04) meshes. The SO3 molecule was placed in 
the middle of a large cubic unit cell (a=20 A) with a single 
k-point. 

In the FPLO approach, we performed band- structure cal- 
culations using the non-relativistic, the scalar-relativistic and 
the full four-component relativistic versions of the full poten- 
tial local orbital minimum-basis band- structure method[23 ] 
(FPLO version 7.00-27). The Local density approximation 



TABLE I. Comparison of the experimental and calculated results for 
the EMF [V] of the lead-battery reaction (1). 



Method 


Level of relativity 


A 


A 




NR SR FR 


FR-NR SR-NR 


BAND VWN +0.55 +2.52 +2.27 
PBEsol-D +0.21 +2.25 +2.02 


+1.72 
+1.81 


+1.97 
+2.04 


FPLO PW92 
PW92b 


+0.41 +2.20 +2.10 
+0.39 +2.21 +2.11 


+1.69 
+1.72 


+1.80 
+1.82 


Av. 
Exp. a 


+0.39 +2.30 +2.13 
+2.107 


+1.74 


+1.91 


a Ref. L24J 



was used l25l for the exchange-correlation functional. Here, 
the optimal volumes of the crystal structures were determined 
by calculating the experimentally determined structures at dif- 
ferent volumes and fitting a 4th-order polynomial equation of 
state to the data points. To sample the first Brillouin zone and 
evaluate the k-space integrals we use a linear numerical in- 
tegration scheme with 24x24x24 (Pb and Sn), 12x12x12 
(PbO, Pb0 2 , SnO, Sn0 2 ) and 6x6x6 (PbS0 4 , SnS0 4 ) k- 
point meshes. The SO3 molecule was calculated without pe- 
riodic boundary conditions. 

The resulting reaction energies for the lead battery reac- 
tion 1), calculated at the FR (fully relativistic), SR (scalar 
relativistic = without spin-orbit coupling), and NR (non- 
relativistic) levels are compared to experiment in Table [I| We 
manage to reproduce the absolute voltage of the lead battery 
reaction within about 0.2V. Taking the four calculations in Ta- 
ble|l]at face value, our calculated absolute voltage for reaction 
(1), in H 2 SO 4 10H 2 O, is +2.13 V while its relativistic part is 
+1.74 V. The relativistic increase of the oxidative power of 
/3-Pb0 2 (s) is indeed the largest contribution. The PbS0 4 (s) 
contribution to EMF follows, and has the same sign. The third 
largest contribution comes from the spin-orbit coupling effects 
in metallic Pb. At 'PBEsol-D' level of theory these three con- 
tributions are +1.58, +0.27 and -0.06 V, respectively, and were 
seen to be essentially method-independent. For details on in- 
dividual contributions, see Fig.[T] 

Comparing Pb with its lighter congener Sn, it has been 
noticed before that no corresponding 'tin battery' exists! 26]. 
This was attributed to the lower oxidative power of tin diox- 
ide, as compared to lead dioxide. We find support for this con- 
clusion in our calculations, since a tin battery would roughly 
correspond to a lead battery without relativistic effects (see 
below). Indeed we found the largest relativistic shifts in the 
lead dioxide. In Fig. [2] we show the FR, SR, and NR den- 
sities of states (DOS) of Pb in /3-Pb0 2 (s). We see that the 
Pb 6s character is evenly distributed between the filled and 
empty states fT6lL making this shell "half-oxidized". Further- 
more, the size of the band gap decreases with increasing level 
of relativity (NR, SR, FR). Note also the relative shifts and 
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changes in the Pb 6s population at different levels of relativ- 
ity. The Pb 6s states are significantly stabilized by inclusion 
of relativity also in /3-PbC>2 . 

We derived above the cell voltage of reaction (1) by using a 
combination of experimental and theoretical results. In order 
to validate the solid-state DFT part - something which cannot 
be taken for granted - we also considered the simplified 'toy 
model' reactions, taking place entirely in the solid-state, 



M(s) + M0 2 (s) -> 2MO(s), 



(5) 



with M being Pb or Sn. The reaction energy of this model 
can be calculated completely by DFT, and should be close 
to experiment. We find that this is the case, with the volt- 
ages for Pb being within 0.05 V of the experimental value for 
the two different DFT implementations and several different 
exchange-correlation functionals, see Fig. [3] This suggests 
that our theoretical calculations accurately describe reality. In 
Fig. [T] we show the calculated energies of formation, Ef, of 
individual species at various levels of relativity. The relativis- 
tic shift is most pronounced for solid /3-PbC>2 . For solid PbO it 
is approximately five times smaller and has the same sign. In 
the case of the tin 'toy-model', the largest relativistic contri- 
bution arises from SnC>2, followed by SnO and Sn, for which 
the relativistic shifts are only slightly different. 

Returning to a qualitative discussion, we mean here by rel- 
ativistic effects anything that depends on the speed of light or, 
more technically, the results of the Dirac versus Schrodinger 
one-electron equation. The main effects are the stabilization 
of all ns and np shells, the destabilization of the nd and nf 
shells and the spin-orbit (SO) splitting of the p,d,f, . . . shells. 
For example, the fully relativistic (FR) and non-relativistic 
(NR) atomic orbital energies of tin and lead are shown in 
Fig. [4] As noticed beforeQ H), the non-relativistic values 
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FIG. 1. The relativistic shifts in energies of formation, Ef (per for- 
mula unit), calculated at DFT/PBEsol-D level. The non-relativistic 
energy, NR = E/(NR), is chosen as reference; Asr = Ef (SR) - 
£/(NR), and A F r = E f (FR) - E f (NR). Note the smallness of the 
effect on SOs(g). 
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FIG. 2. Total and partial DOS (states/eV/cell) of Pb in /3-Pb0 2 cal- 
culated at DFT/VWN level (the Fermi energy, Ef, is set to 0). Note 
the pronounced relativistic shift in both occupied and unoccupied Pb 
6s-states. (NR - non-relativistic, SR - scalar-relativistic, FR - fully- 
relativistic). Orange: Pb 6s, green: Pb 6p. 
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FIG. 3. Effect of relativity on the EMF of model-reaction, M(s) + 
M0 2 (s) -> 2MO(s), calculated with four different methods (A-D). 
Experimental values for M=Pb (+0.82 V) and M=Sn (-0.085 V) are 
indicated with the horizontal dashed and dotted lines, respectively. 
A) BAND (VWN), B) BAND (PBEsol), C) BAND (PBEsol-D), and 
D) FPLO (PW92). The level C has dispersion corrections. 



are similar for Sn and Pb while the relativistic ones are not. 
Qualitatively, the tendency of Pb to be predominantly diva- 
lent can be related to the increased binding energy of the Pb 
6s shell. Relativistic effects have also been suggested as re- 
sponsible for the changing the structure of metallic lead from 
diamond to fccl9ll27lL as well being necessary to determine 
phase transitions between fee, hep, and bec structures |[T0ll . 

Moreover, it has been shown that the valence- shell relativis- 
tic effects of the Periodic Table scale roughly as Z 2 , Z being 
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FIG. 4. Relativistic (FR) and non-relativistic (NR) Hartree-Fock or- 
bital energies for the tin and lead atoms in their ns 2 np 2 ground state. 
Note the relativistic stabilization of the ns level (lowest in the figure), 
leading to higher oxidative power for Pb(IV) than for Sn(IV). Data 
from Desclaux|28 1. 



the full nuclear charge 0. For tin and lead the ratio is 



'Z(Sn)" 


2 


"50" 


[z(Ph)\ 




_82_ 



0.372, 



(6) 



while the fully relativistic voltage changes, AE, calculated 
for the 'toy models' ([5]) yield 



' AE(Sn) ' 


2 


"0.34V 


_A£(Pb)_ 




0.86V_ 



= 0.395. 



(7) 



Indeed, the EMF calculated for the non-relativistic lead 'toy- 
model' is similar to the analogous tin model reaction, treated 
at fully-relativistic level. It also explains why the analogous 
reaction (1) for tin is unknown. 

Concluding, the lead-acid battery belongs to those famil- 
iar phenomena, whose characteristic features are due to the 
relativistic dynamics of fast electrons when they move near a 
heavy nucleus. In this case the main actors are the 6s valence 
electrons of lead, in the substances involved. This insight may 
not help one to improve the lead battery, but it could be useful 
in exploring alternatives. Finally, we note that cars start due 
to relativity. 
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